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ABSTRACT 

The effect of mobility ratio on areal sweep efficiency and con- 
ductance ratio has been made for direct and staggered line drive well 
patterns* This study was made using miscible fluids, artificially 
consolidated sandstone models, and a photographic method of data 
recording. The mobility ratio was varied between 0. 1 and 12, which 
is the range common to water flood operations. 

A method of predicting the performance of waterfloods in 
stratified reservoirs using conductance ratios is treated. This 
technique is not limited in number of strata and indicates one method 
of considering a free gas saturation. 
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INTRODUCTION 



The importance of secondary recovery techniques in petroleum 
engineering is continually increasing and is taking an increasing per- 
centage of reservoir research activity. Since primary reservoir deple- 
tion normally will leave approximately 75 per cent of the original oil 
in place, the reason for this increasing role of secondary recovery is 
obvious . 

Perhaps the oldest and most used method of recovering oil by 
secondary techniques is the waterflood. This involves injection of 
water into the reservoir to displace the oil. Water injection requires 
expenditure of energy and funds, therefore, reliable methods for pre- 
dicting waterflood results are necessary. 

Extensive study has gone into conformance factor or sweep 
efficiency predictions for various well patterns at varying mobility 
ratios. This permits predictions for single stratum reservoirs. 

Most reservoirs are composed of several strata, however, and 
performance predictions for stratified reservoirs are not feasible 
with only sweep efficiency data available. In predictions for stratified 
systems, it is necessary to place performance on a time basis. This 
requires that flow potential or flow rate for a pressure drop be avail- 
able for the system under investigation. 
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Flow potential data (conductance ratios) are available for all 
7 , 20, 34 

regular patterns with the exception of the line drive. The pur- 

pose of this study was to obtain line drive conductance ratios and to test 
methods for approximating conductance ratios analytically. A method 
of predicting performance of a stratified system using these data is 
also presented. 



WATER FLOOD PREDICTIONS 



In studies of petroleum reservoirs, much qualitative and quanti- 
tative use has been made of models to more completely understand the 
phenomenon of oil production. Since no two reservoirs are exactly the 
same, model studies are limited in scope. Therefore, in a rigorous 
interpretation, models can only be used for studies of a particular 
field; however, by generalizing parameters, study of a reservoir seg- 
ment will yield data for qualitative application in a number of situations. 

A great deal of work has gone into the determination of scaling 
criteria for model studies. Probably the most 

complete analysis covering general situations has been presented by 

17 

Geertsma, Croes, and Schwarz. A strict interpretation of this 

analysis, however, leads one to the conclusion that the model must 

be identical to the prototype in every respect. Since this is completely 

impractical, it is necessary to decide which reservoir factors are of 

greatest importance and construct the model accordingly. A detailed 

description and application of model scaling for edgewater drive reser- 

4 

voirs has been presented by Caudle. 

The efficiency of a secondary recovery project, in which a fluid 
(usually water) is injected to displace the oil in the reservoir, can be 
defined by the volumetric sweep efficiency. In reservoir studies, 
models are normally two dimensional representations and the displace- 
ment efficiency may be defined in terms of areal sweep efficiency. This 
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is the ratio of the area swept by the injected fluid to the total model area. 

Generally, the areal sweep efficiency is determined by the well 
pattern geometry and the mobility ratio. This assumes that the reser- 
voir is homogeneous and that gravitational forces are not present. The 
mobility ratio is defined as the ratio of the sum of the mobilities of the 
fluids behind the front to the sum of the mobilities of the fluids ahead 
of the front. The mobility {\) of a fluid in a porous media is expressed 
as the effective rock permeability (k) to the fluid of interest at its 
saturation divided by the fluid viscosity (\l ) . 



(k/|u) swept 
(k/ \j) unswept 




Mobility Ratio 



( 1 ) 



Early model studies utilized the electrolytic model to represent 

waterfloods. This model uses the principle of ionic flow to simulate 

fluid flow. Use of this model is limited to mobility ratios of one. 

27 

Muskat, Wyckoff, and Botset were among the first to utilize this 

method to determine breakthrough sweep efficiencies for various well 

26 

patterns. In a later publication, Muskat and Wyckoff analyzed the 

27 

previous electrolytic model results and concluded that the shape of 
the flood front and the breakthrough areal sweep efficiency were depend- 
ent on the geometry of the system and not on pressure differential, well 

24 25 26 

spacing, or rock properties. It is noted that Muskat ’ * indicates 
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that the results obtained in the electrolytic models do not consider the 
effects of two fluids with different viscosities. In other words, the elec 
trolytic model is valid for a mobility ratio of one only. 

The electrolytic model and the similar but more sophisticated 

potentiometric model have since been used in the study of various reser 

voir conditions. Limited use has also been made of heat flow models to 

22 

simulate fluid flow with a mobility ratio of one. 

Sweep Efficiency 

Early studies on the effect of mobility ratio on areal sweep effi- 
ciency utilized potentiometric models. Aronofsky"'’ reported results of 
a potentiometric model study of mobility ratio effects. A stepwise pro- 
cedure was utilized to vary the ratio of the depth of the electrolyte 
behind the front to that ahead of the front. This variation in depth made 
the electrical conductivity ahead of the front different from that behind 
the front, thereby simulating fluids of different mobilities flowing in a 
porous media. After computing a new frontal position, the electrolyte 
depth was changed to fit the new frontal position and the process was 
repeated. In this iterative manner, sweep efficiencies for mobility 
ratios other than one were^approximated by advancing the front in a 
series of steady-state steps. 

The effect of mobility ratio on areal sweep efficiency has also 

28 



been studied with an electric analog model by Nobles and Janzen. 
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This study utilized a resistance network and an iterative process similar 
to Aronofsky* s* to simulate frontal movement. 

g 

Cheek and Menzie studied the effects of mobility ratio on areal 
sweep through the use of a fluid mapper. This model uses fluid flow 
between two closely spaced parallel plates to simulate fluid flow in a 
porous media. Results are reported for five spot and direct line drive 
well patterns. 

32 

Slobod and Caudle introduced the X-ray shadowgraph technique 
to study fluid flow phenomenon. This new technique used a porous plate 
and fluid flow with one of the fluids containing an X-ray absorbing mate- 
rial so that a contrast between displaced and displacing fluid could be 
observed on X-ray shadowgraphs taken during the experiment. By 
taking successive X-ray exposures, a complete history of the flood was 
obtained. 

13 

Dyes, Caudle, and Erickson used the X-ray shadowgraph 

method to illustrate that a significant amount of production could take 

place after initial breakthrough. The results of this study demonstrated 

that, with most water floods, an ultimate areal sweep efficiency of 85- 

100 per cent could be obtained for various mobility ratios at an economic 

water-oil ratio. The only previous reports which indicated consideration 

21 

of production after breakthrough were made by Hurst for a mobility 

1 5 

ratio of one and Fay and Prats for a mobility ratio of 4. 0 using 



numerical techniques. 
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The X-ray shadowgraph has been used extensively to study a 

5 , 6 , 7 , 11 , 13 , 14, 16 

variety of reservoir conditions at varying mobility ratios. 

Included among these studies are five -spot pattern floods, stratified 

five-spot floods, and the effect of fractures on flood performance. 

1 8 

Habermann utilized a method similar to the X-ray shadow- 
graph for studying fluid flow. This method consisted of artificially 
consolidated sand models covered with clear plastic. One of the fluids, 
either displaced or displacing, contained a dye which permitted visual 
observation of the front as it progressed. Photographs taken during 
the flood permitted calculation of areal sweep efficiency. 



Conductance Ratio 

The use of areal sweep efficiency data permits ultimate re- 
covery predictions for single stratum reservoirs. In order to predict 
performance of stratified reservoirs and also to predict recovery of 
single stratum reservoirs on a time basis, it is necessary to know the 
change in the flow potential or flow rate per unit pressure drop for the 
system under study. 

24 

Prediction of flow potential was first discussed by Muskat in 
terms of conductivity (production rate per well per unit -pres sure 
differential). This permitted predictions for systems with mobility 
ratios of one. 

2 

Aronofsky and Ramey first introduced the concept of conductance 
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ratio in a potentiometric study of mobility ratio effects on a five -spot 
water flood. The conductance ratio (y) is defined as the flow rate (q) 
per unit pressure drop (AP) between the injection and production wells 
at a specified time to that conductivity with only the original fluid pre- 
sent. 



= Conductance Ratio (2) 

Analytical methods of studying the change in flow potential has 

been presented in terms of dimensionless injectivity by various methods. 

30 

Prats, Matthews, Jewett, and Baker have computed injectivity for 
five -spot systems by considering that the reservoir consists of several 

19 

layers, each having a uniform but different permeability. Hauber 
presented another approximate method which considers that flow uti- 
lizes the stream tubes and iso -potential lines obtained by the solution 

of differential equations for a mobility ratio of one. 

12 

Deppe presented a method to approximate conductivity for vary- 
ing mobility ratios and for varying well patterns. In the five -spot 

12 

pattern used as an example, Deppe essentially divided the pattern 
into two regions of radial flow. The first region was around the injec- 
tion well and covered almost one -half of the area. The second radial 



(q/AP) time t 
(q/AP) initial 



flow region was the same except that it was around the production well. 
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The results of Deppe's study provided injectivity (conductance 
ratio) curves which were not too unreasonable for five -spot and nine- 
spot well patterns. Closer agreement was obtained for unfavorable 
mobility ratios, which is normally true when using this type of approxi- 
mation methods. 

In an attempt to improve on the work of Deppe, this study has 
divided the pattern into 3 regions as illustrated by Figure 1 for a 
staggered line drive pattern. The first region surrounds the injection 
well where the largest pressure drop is encountered. This region is 
considered to have radial flow and has a radius equal to one -fourth the 
well row spacing. The second region is considered to be linear, with 
only slight pressure change, and extends to the areal sweep at break- 
through. The third region is one of radial flow around the production 
well where the pressure drop is considered to be a function of frac- 
tional flow from the swept region. 

Conductance ratio for the system is computed by considering 
a constant flow rate which reduces the equation to 



Y 




(3) 



Pressure drop across the system is obtained by summing the pressure 
drops in the three regions at the time of interest. 



10 




Injection Well 



ri Radius of radial flow area surrounding injection well 

rp Radius of radial flow area surrounding producing well 

0C Angle open to flow from swept region after breakthrough 



FIGURE 1 



STAGGERED LINE DRIVE -- ASSUMED FLOW REGIONS 
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With only single fluid flow in the system (or at initial conditions) 
the pressure drop in each region can be represented by modifications of 
Darcy's Law as 



AP 



Ai 



q In ri/ rw 

X T.'O'tfh""' 



(4) 



AP 



Bi 



q L/w 



X 1. 127h 



(5) 



AP 



Ci 



q In rp/ rw 

~irr7mr 



( 6 ) 



which when added, gives an initial pressure drop for the system of 



A Pi = ^ 7^ frgh U n r i/ rw + 6.28 L/w + In rp/rw) (7) 
u 

After the flood has started and two fluids are flowing in the model, the 
pressure drop can be represented as a summation of the pressure drop 
in each fluid region. When the front is in the first region the pressure 
drop therein can be represented as 

q In rf/ rw , q In ri/rf 

AP At X. 7. 08h + X 7. 08h 

s u 

q (1/M) In rf/ rw + In ri/rf 

~ X 7. 08h 

u 



( 3 ) 
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When the front is in the second region the pressure drop therein is 



After breakthrough, the flow rate in the third region must be summed 

for parallel flow since two fluids are flowing into the producing well. 

Consider that portion of the well illustrated in Figure 1 or one-quarter 

of the well. Let a equal the angle open to flow from the swept region 

and (rr/2-a) equal the angle open to flow from the unswept region. 

30 

Then the fractional flow from the swept region is 




qi , q(L-i) 



\ 1. 127hw + \ 1. 1 27hw 



s u 




( 9 ) 



u 



X a 
s 



(10) 






then 



a 



fs rr/2 



(ii) 



M-fs(M-l) 
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The pressure drop in the swept and unswept region is 



AP 



Cs 



fs q In rp/rw tf/2 

7.08'K \"'< i 

s 



( 12 ) 



AP 



Cu 



(1-fs) q In rp/rw tt/2 

7. OSh X (rr / 2- a) 
u v 



(13) 



For parallel flow, the pressure drop given by equation (12) must equal 
that of equation (13). Therefore, the pressure drop in the third region 
may be obtained from either the equation for the swept or the unswept 
portion. By substituting equation (11) into (13) we obtain 

(1-fs) q In rp/ rw 

AP Ct = 7.08h \ u [l-fs/[M-fs(M-l)]j (14 



In summing equations (8), (9), and (14) the following conditions 

apply: 

1. When Es<Es^: jj/w = 0, fs = 0, and rf = f(Es). 

2. When Es A <Es<ESg^,: rf = ri, fs = 0, and j£/w = f(Es). 

: rf = ri, A/w = L/w, and fs = f(Es). 



3. When Es>Es 
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The total pressure drop at any time is then 



AP t = ^ y Q8h In rf/rw + In ri/rf 

u * \ 

+ 6. 28 jjji- - l^j X/w + 6. 28 L/w 



(15) 



(1-fs) In rp/rw 



The conductance ratio from equation (3) becomes the ratio of equation 
(7) to equation (15) or 



y = 6P^/ APj. = (In ri/ rw + 6. 28 L/w + In rp/ rw)/ 



(1/M) In rf/rw + In ri/rf + 6. 28 






4/w 



(16) 



, L « T , (1-fs) In rp/rw \ 

+ 6.28 L/w + l_ f s/[ M _fs(M-l)V 



Values for L/w were obtained by using both arbitrary values and by 
assuming that the pressure drop across a linear system is equal to 
that across a line drive with well radius equal to ri. Muskat's equation 
for flow rate in a line drive is: 



q(i d/a - 1. 17 + 0. 637 log a/ri 

kh 2. 254 



AP 



(17) 
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Setting equation (17) equal to (5) and solving for L/w yields 



L/w 



d/a-1 . 17+0.637 log a/ ri 

176 



(18) 



which is valid when d/a>l. Calculations were performed by programming 
the above for the Controlled Data Corporation 1604 digital computer. 
Appendix D contains the flow diagram and FORTRAN language program 
utilized in these computations. A conversion table for the FORTRAN 
symbols is also provided in Appendix D. 

In the calculations for conductance ratios, various values of L/w 

were assumed in an attempt to match the model study results. Using 

equation (18), L/w was 0.357. Other values used for L/w were 0 and 

0. 1. As can be seen from the plots of Appendix D, varying L/w made 

only small changes in the conductance ratio. Comparison plots for L/w 

of 0. 1 and as obtained from the model studies, indicate results similar 

12 

to those obtained by Deppe. Agreement is fair for unfavorable mobility 
ratios and not too good for favorable mobility ratios. 

This method can be used for rough approximation of conductance 
ratios for varying well patterns by modifying the initial values in the 
program to correspond to those of the pattern under study. However, 
it is believed that fluid flow model study results would be more reliable. 

The effect of mobility ratio on conductance ratio has been 
studied for various well patterns using fluid flow models. Caudle and 
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Witte have published data for the five -spot, the inverted skewed seven- 

20 

spot has been studied by Hickman, and the inverted nine -spot was 

34 

studied by Watson. This study provides data for line drive well 
patterns . 

Stratified Reservoirs 

In a reservoir composed of strata of varying permeabilities 
and porosities, performance prediction becomes more complex than 
for a single stratum. To illustrate, consider a two strata system of 
equal characteristics with fluids of equal viscosity and with the perm- 
eability in the first stratum twice that of the second. From Darcy's 
Law 



k AP 
U 



(19) 



it can readily be seen that the first stratum will have a frontal advance 
twice that of the second strata and that at breakthrough for the first 
stratum the second stratum will have a volumetric sweep of only half 
of the breakthrough sweep efficiency. In this simple example, a mis- 
cible system was assumed in order that mobility ratios of one could 
be used for each strata. For mobility ratios other than one, the linear 
relationship of flow velocity to permeability no longer exists in a 
stratified system and then becomes a function of distance traveled or 



area swept. 
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In order to apply the data obtained in this study to a stratified 
reservoir, it is necessary to put the performance of each strata on an 
equal time basis. This is most easily done through the use of the con- 
ductance ratios obtained in this and other studies. It is noted that this 
method is developed and accurate for separated strata or for strata 
without cross flow. When cross flow between strata exists and the 
mobility ratio is greater than one, this method provides a good 
approximation for predicting water flood results. A similar method 

using theoretical dimensionless injectivities has been presented by 
30 

Pratts, et al. The method developed herein permits performance 

prediction as accurately but more simply than the method described 
30 

by Prats , et al. 

To develop this time relationship, consider Darcy’s Law (19) 

and let 

v = q/A 

and 

A = Vd/i 

Darcy's Law can then be written in terms of conductance as 

q/AP = — ir-TT = conductance (20) 

U l M 

Consider now two adjacent parallel strata in a reservoir with single 
fluid flowing; the variables between these two strata are permeability (k) 
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and displaceable volume (Vd). 

Therefore the conductance between a stratum 1 and any stratum 
j will vary as their permeability and displaceable volume ratios vary 
or 

q. k.Vd. q, 

J ^ J J -1 / 9 i 

apt - Tqw L "ap^ (21 

By definition, 



Y = (q/AP) t / (q/AP). 



( 22 ) 



then 




( 23 ) 



and 




( 24 ) 



solving for dt, integrating, and letting c equal 




ct 




dWi 



( 25 ) 



19 



Other equations needed to express stratified reservoir perform- 
ance are for total flow from the swept region (fs) and total conductance 
ratio (y). 

By definition, 

fg = qs/q = E qs. / E (26) 

then 



and 



from (7) 



qs . = fs. q. 

J J 



q j = Y j Cq A ZSF7T, 



k.Vd. 

<Vi = C q l\ 



then by substituting into (26) we have 

E v. fs. k. Vd. 

fs = J J J J 

1 " E y. k. Vd. 



(27) 



(28) 



(29) 



(30) 



'J J J 

The expression for overall conductance ratio can be derived in this 



same manner as 



T V- k. Vd. 
_ Y J J J 
Y - 2 k. Vd. 

J J 



(31) 



An example performance prediction for a stratified reservoir 
using this method is described in Appendix E. A possible method of 
considering an initial free gas saturation is also mentioned in Appendix 
E. 



EXPERIMENTAL APPARATUS 



Model Theory 

The majority of oil reservoirs are drilled with a regular pattern 
of wells. Probably the most common pattern utilized in this country is 
a square array of wells on a 40 acre spacing. This well pattern, illus- 
trated by Figure 2, readily lends itself to secondary recovery operations 
normally referred to as the Direct Line Drive, which is obtained by con- 
verting alternate rows of wells to injection wells. 

The Staggered Line Drive is also obtained by converting alternate 
rows of wells to injection wells; however, the wells are not drilled in a 
square array but rather in an offset fashion. The well pattern for the 
Staggered Line Drive system is shown by Figure 3. 

With the regular well systems shown by Figures 2 and 3, it can 

readily be seen that there will exist numerous mirrow images divided 

by lines of symmetry. These lines or planes of symmetry represent 

invarient streamlines or boundaries across which no flow will occur, 

9 

as long as production and injection quantities are uniform. Since the 
planes of symmetry are in effect flow boundaries in a homogeneous 
system, it is apparent that only the smallest segment bounded by 
symmetry planes need be considered. Figures 4 and 5 show enlarged 
views of the model and its relationship to the line drive systems studied. 
As can be seen, the model utilized represents a line drive pattern where - 
in the row spacing equals the well spacing. 



20 
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O Production Well 

9 Injection Well 

Lines of Symmetry 



FIGURE 2. DIRECT LINE DRIVE WELL ARRAY 
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FIGURE 3. STAGGERED LINE DRIVE WELL ARRAY 
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FIGURE 4 . DIRECT LINE DRIVE PATTERN 
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FIGURE 5. STAGGERED LINE DRIVE PATTERN 
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Model Construction 



The model utilized in this study was an artificially consolidated 

sand plate constructed generally as described by Watson*^ and Hickman. ^ 

The artificial sandstone was composed of a clear epoxy resin and 

60-100 mesh Ottawa sand. Actual construction of the artificial sandstone 

34 20 4 

has been reported in detail by Watson, Hickman, and Caudle. In 
order to avoid any possible loose streaks which may have occurred during 
packing of the artificial sand a minimum of 2 inches was cut from the 
model sides by using a diamond saw. To seal the model, an extremely 
viscous epoxy resin was spread over all exposed sand faces. The high 
viscosity of the epoxy resin was necessary to prevent imbibition of the 
coating into the model and was obtained by applying the coating just 
prior to it being completely n set-up. M 

Injection and producing '’wells' 1 were constructed from 1/2 inch 
thick acrylic plastic plate as shown by Figure 6. Swagelok fittings for 
1/8 inch tubing were utilized to connect the "wells" to the injection 
systems. Prior to attaching the plastic plates to the model, well bores 
were reamed out at the well corners. 

As seen from Figure 6, the model utilized was 32 cm in length, 



16 cm in width, 0. 6 cm thick, and had a well bore diameter of 0. 2 cm. 
This represents a d/rw ratio of 320. 
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32 cm 





FIGURE 6. LINE DRIVE MODEL 
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Fluids 

The fluids utilized in this study were miscible in all proportions. 
The use of miscible fluids to simulate reservoir fluids in model studies 
of mobility ratio effects has been well established with experimental data 
by Slobod and Caudle‘S and in other publications. ^ 9 ^ 9 ^ 9 ^ 

By using miscible fluids in this study, there are no effects of 

capillary forces as will be found in immiscible systems. The absence 

of capillary effects eliminates the problem of saturation gradients in 

determining the effective mobility ratio. Also by utilizing miscible 

fluids, the effective permeability to each fluid is the same, thereby 

eliminating the need to determine relative permeabilities. This system 

therefore models a water flood with little or no flow of oil behind the 

swept — then reduces to the 
unswept 

viscosity ratio of U1 ? . SW ^P* Therefore, only the viscosity of the 
7 ((i) swept 7 7 

fluids needed to be changed in order to vary the mobility ratio of the 
system. 

To obtain fluids of the desired viscosities, a mixture of gly- 
cerine and water was used. The quantity of glycerine added for the 
desired viscosity was determined on a weight per cent basis through 
the use of standard chemical tables and checked with an Ostwald 
viscosimeter. 

In the photographic technique used in this study, one of the fluids 



water front. The mobility ratio, 



(k/u) 

(k/|i) 
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contained a dye to establish a contrast between the displacing and the 
displaced fluids. Since the epoxy used to consolidate the sand renders 
the grains oil wet, oil soluble dyes could not be used. It was found that 
an ordinary household food coloring used in about 5 per cent volumetric 
proportions with the water obtained the desired contrast. 

Injection System 

20 

The constant rate pumps, utilized and described by Hickman 
34 

and Watson were used in this study. This pump is basically a hydraulic 
cylinder using a screw type mechanism to move the piston. It is driven 
by a gear and chain arrangement from an electric motor through a 
variable speed reduction gear. By changing the gears in the drive mech- 
anism and/or the speed output of the reduction gear, the desired injection 
rate could be obtained. 

A pressure gauge was located at the injection well. The producing 
well was permitted to discharge at atmospheric pressure, thereby per- 
mitting the gauge to indicate the pressure drop across the model. 

Data Recording 

A 35 mm camera was used to record the frontal position, pressure 
reading, and time during the runs. The camera was actuated by a D. C. 
timing motor which advanced the film and made one exposure per minute. 
Exposures could be made at longer intervals by using a "sign flasher" to 
interrupt current flow to the timing motor. This was done during the 
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slower runs. A red filter was used to obtain better contrast between 
the blue and clear fluids on black and white medium speed panchro- 
matic film. A more complete description of the photographic and 

4 

developing technique has been previously written by Caudle. 

Figure 7 is a schematic diagram of the equipment used in this 

study. 

Typical photographs are shown in Figures 8 and 9. Figure 8 
is an exposure of a direct line drive at a mobility ratio of one and Figure 
9 is for a staggered line drive at a mobility ratio of twelve. Both expo- 
sures indicate frontal positions shortly before breakthrough. 



Production 
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FIGURE 7. SCHEMATIC DIAGRAM OF EQUIPMENT 
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FIGURE 9. PHOTOGRAPH OF MODEL DURING A STAGGERED 
LINE DRIVE RUN, M = 12. 







EXPERIMENTAL PROCEDURE 



The performance of both direct and staggered line drive patterns 
were studied for mobility ratios between 0. land 12 to determine the 
effect on areal sweep efficiency, fractional flow, and conductance ratio. 

The model was initially evacuated by means of a vacuum pump. 
After complete evacuation, the model was saturated with water. At 
this time, all lines and the pressure gauge were connected to the model 
and the injection pump (Figure 7). The injection pump was loaded with 
dyed water as the initial displacing fluid, the model, pressure gauge, 
and clock were placed on a light table beneath the 35 mm camera and 
the direct line drive runs were started. 

As injected fluid moved into the model, the automatic camera 
recorded the advancing front positions and injection pressures at regular 
intervals. The camera was stopped after approximately 95 per cent of 
the model had been swept. In either case, injection was continued until 
complete displacement of the "oil analog" had occurred. 

In the second run, the injection system was loaded with fluid of 
higher viscosity and the above process was repeated with the clear 
viscous fluid displacing the dyed fluid. By alternating injection between 
clear fluid of varying viscosity and the dyed fluid, the range of mobility 
ratios from 0. 1 to 12.0 was covered in five runs. 

After completion of a set of runs covering the range of mobility 
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ratios under study, the roles of the injection and producing wells were 
reversed and the five runs were repeated. This procedure was then 
followed for the staggered line drive system. The two runs at the 
same conditions were made in opposite directions to average out any 
difference in the pressure drop around the individual wells. 

The flow rate was set prior to the start of each run and was 
maintained constant throughout the run. The rate was chosen so that 
the injection pressure would not exceed the gauge scale and so that a 
sufficient number of photographic exposures would be made during the 
run to permit interpretation. The areal sweep efficiency is independent 
of flow rate in the line drive floods. Since the conductance ratio 
(Equation 2) is a dimensionless ratio, the actual rate utilized is of no 
concern. It is important, however, to ensure that the flow rate remains 
constant during each run so that the time scale can later be converted to 
pore volumes of water injected and the conductance ratio computed on 



the ratio of pressure drop. 



ANALYSIS OF DATA 



Areal Sweep Efficiency 

All data was recorded on film during each run. After developing 
the film, a microfilm reader was used to enlarge the recorded image. 
Tracings were then made of the model and the frontal positions during 
sweepout. Appendix A contains typical tracings of frontal positions at 
various mobility ratios for both the direct line drive and the staggered 
line drive. The area of the model and the area behind the frontal posi- 
tion were obtained with a polar planimeter. The areal sweep efficiency 
was then obtained for each frontal position by dividing the area swept by 
the area of the model. Tables of data are contained in Appendix B. 

Areal sweep efficiencies were plotted versus an elapsed time 
scale. Since runs at the same mobility ratio, but in opposite directions, 
were not necessarily made at the same rate, a time correction factor 
was obtained by plotting areal sweep data for the second run on the 
first curve obtained and then noting the time correction necessary which 
would make the curves most closely coincide. This permitted easy 
correlation of data for runs of the same mobility ratio made in opposite 
directions . 

Since each run was made at a constant rate, it is easily seen that 
the time scale of the basic data plots is directly related to the volume of 
displacing fluid injected by a constant ratio for each particular rate. 
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Prior to breakthrough of injected fluid, a plot of areal sweep versus 
volume injected is a straight line. By noting the area swept and the 
elapsed time at breakthrough, all subsequent time values may be con- 
verted to displaceable volumes of injected fluid by the ratio of area 
swept to elapsed time at breakthrough. Appendix B contains plots of 
all mobility ratios studied for the two well patterns considered, with 
tables of data derived from these plots. 

Plots of areal sweep versus displaceable volumes injected can 

be found in Appendix C. The displaceable volume injected is a more 

useful term than actual volume (or time at constant injection rate) due 

13 

to its dimensionless nature. Dyes, et al has defined displaceable 
volume as the product of the hydrocarbon pore volume of the unit 
pattern and the volumetric displacement efficiency. 

It is noted that areal sweep data from the first 5 runs did not 
correspond to other runs for the direct line drive made in the opposite 
direction as regards areal sweep efficiency. After study of the photo- 
graphed flood front shapes, it was noted that an air pocket had been 
present near the production well during these runs, thereby distorting 
streamlines and causing invalid areal sweep data. This error was 
noted after completion of all runs and an attempt was then made to re- 
make these runs. However, the injection well became "plugged" during 
these "reruns" and areal sweep data for the direct line drive for mobility 
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ratios of 3. 13, 0. 106, and 12. 1 is based on one run only. The cause 
of this "well plugging" is attributed to impurities in the injected fluid, 
which could have been prevented by installing a filter on the injection 
line. This air pocket did not interfere with runs subsequent to run 5 
since a number of days elapsed between run 5 and 6 and the model was 
routinely evacuated before making run 6. 

Conductance Ratio 

2 

The conductance ratio is defined by Aronofsky and Ramey 
(q/AP) / (q/AP). which reduces to Ap. / £P when the flow rate is 
maintained constant. Pressure data taken from the injection well pro- 
vides AP at a particular areal sweep efficiency and permits direct 
calculation of the conductance ratio. 

From the above equation, it can be seen that if the pressure 
increases during a run, the conductance ratio will decrease and a de- 
creasing pressure will cause an increasing conductance ratio. The 
conductance ratio decreases if the mobility ratio is favorable and in- 
creases if the mobility ratio is unfavorable. At a mobility ratio of one, 
the conductance ratio will remain constant during a flood. 

Plots of conductance ratio versus areal sweep can be found in 
Appendix C. 

The ratio of the distance between injection and production wells 

/ 2 
to the wellbore radius d/ rw has been shown by Aronofsky and Ramey 
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to effect the conductance ratio. Although the model was constructed with 

a d/ rw ratio of 320, slight differences in the well bore radius and plugging 

around the well bore will tend to change this ratio. 

The steady state flow equation for single phase flow in line drive 

24 12 

patterns was derived by Muskat and presented by Deppe in practical 
units as 



0.001538 k/iih AP 

0. 682 d/a -0. 798 + In d/ rw 



(32) 



This equation is valid for both direct and staggered line drive well 
patterns with d/a>l. 

Using the information obtained from the model runs, the effec- 
tive d/rw ratio for the direct line drive was found to be 922. The 
effective d/rw ratio for the staggered line drive was calculated as 1110. 
These values for d/rw are somewhat smaller than those normally en- 
countered in the field. However, it is considered that reasonable 
accuracy will be obtained by using the conductance ratios obtained in this 
study without correction. 

If it is desired to correct the conductance ratios obtained in this 

study to the correct value for other d/rw ratios, using the method pre- 

2 

sented by Aronofsky and Ramey for line drives yields the following re- 



lationship: 
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Yf 



0, 682(d/a)^- 0. 798+ln(d/ rw ) m + 0 • 434 In ^ rw ^ m 
1/V m [0. 682(d/a) f -0. 798 + ln(d/rw) m ] 



(33) 



where the subscripts f and m denote field and model respectively. 
Fractional Flow 

Up to the time of breakthrough, the fractional flow from the 
unswept region (oil cut) is one. Also up to the time of breakthrough, 
a plot of areal sweep versus displaceable volumes injected yields a 
straight line at a slope of one. It can also be seen that a plot of areal 
sweep versus elapsed time, with injection at a constant rate, will 
yield a straight line prior to breakthrough. After breakthrough, the 
areal sweep curve deviates from a straight line. 

Previous publications^’ ^ have shown that the fractional flow 
from the producing well can be obtained from the slope of the areal 
sweep efficiency curve versus displaceable volumes injected. Since 
elapsed time and displaceable volumes are directly proportional, the 
fractional flow from the unswept region can be obtained by the ratio 
of the slope of the areal sweep efficiency curve at the time of interest 
to the slope of the straight line portion. Tables of fractional flow are 
given with the basic data plots in Appendix B and are plotted versus 
areal sweep efficiency in Appendix C. 



INTERPRETATION OF RESULTS 



Tables of basic data from the model runs are presented in 
Appendix B. Appendix B also includes curves of areal sweep efficiency 
versus elapsed time. Each figure also includes tables of displaceable 
volumes injected and fractional flow data as derived from these curves. 
Plots of displaceable volumes injected, fractional flow, and conductance 
ratio versus areal sweep efficiency for the mobility ratios studied are 
provided in Appendix C. 

Results of this study were compared with those obtained by the 

1 3 

X-ray shadowgraph technique and were found to have a maximum 

deviation of 7% and an average difference of less than 2%. Comparison 

of areal sweep efficiencies at breakthrough and at 95% fractional flow 

13 

from the swept region as obtained by Dyes, et al and as obtained in 
this study are listed in Table I for the direct line drive and in Table II 
for the staggered line drive. 

In order to make the results of this study useful for mobility 
ratios other than those specifically used here, Figures 10 through 15 
indicate the effect of mobility ratio on performance. 

In comparing the results obtained for the direct line drive and 
the staggered line drive. Figure 16 shows that areal sweep efficiencies 
are significantly higher for the staggered line drive at breakthrough. 
However, at a fractional flow from the swept region of 95 per cent 
(assumed economic limit) the areal sweep efficiencies are similar for 
the two patterns at abandonment up to a mobility ratio of approximately 
10. At mobility ratios above 10, previous studies have shown that the 
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41 



areal sweep efficiencies decrease rapidly in this region for a given 
fractional flow. 
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TABLE I 

COMPARISON OF AREAL SWEEP EFFICIENCIES 
FOR DIRECT LINE DRIVE 





Breakthrough 
Sweep Efficiency 


95% Fractional Water Flow 
Sweep Efficiency 


Mobility 

Ratio 


__ „ ,13 

Dyes, et al 


This Study 


_ . ,13 

Dyes, et al 


This Study 


0. 106 


0.92 


0. 85 


0.995 


0.98 


0. 274 


0. 76 


0. 76 


0.99 


0.97 


1.0 


0.55 


0. 55 


0.98 


0.96 


3. 13 


0.46 


0. 46 


0. 97 


0.95 


12. 1 


* 


0. 41 


* 


0. 77 



*Not available 
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TABLE II 

COMPARISON OF AREAL SWEEP EFFICIENCIES 
FOR STAGGERED LINE DRIVE 



Mobility 

Ratio 


Breakthrough 
Sweep Efficiency 


95% Fractional Water Flow 
Sweep Efficiency 


1 3 

Dyes, et al This Study 


Dyes , et al^ ^ 


This Study 


0. 106 


0.91 


0.89 


0.995 


0.99 


0. 274 


0. 85 


0.87 


0.99 


0. 98 


1.0 


0. 75 


0. 77 


0.98 


0.98 


3. 13 


0. 61 


0.66 


0. 95 


0. 96 


12. 1 


* 


0. 55 


❖ 


0. 89 



*Not available 
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FIGURE 10. DIRECT LINE DRIVE -- VARYING DISPLACEABLE VOLUMES for AREAL 

SWEEP versus MOBILITY RATIO 
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FIGURE 11. DIRECT LINE DRIVE -- VARYING FRACTIONAL FLOWS for AREAL SWEEP 

versus MOBILITY RATIO 
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FIGURE 12. DIRECT LINE DRIVE -- VARYING AREAL SWEEP for CONDUCTANCE 

RATIO versus MOBILITY RATIO 
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FIGURE 13. STAGGERED LINE DRIVE -- VARYING DISPLACEABLE VOLUMES 

for AREAL SWEEP versus MOBILITY RATIO 
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FIGURE 14. STAGGERED LINE DRIVE -- VARYING FRACTIONAL FLOWS for AREAL 

SWEEP versus MOBILITY RATIO 
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FIGURE 15. STAGGERED LINE DRIVE -- AREAL SWEEPS for CONDUCTANCE RATIO 

versus MOBILITY RATIO 
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CONCLUSIONS 



The results of this study of line drive flood patterns indicate 
that mobility ratio has a definite effect on the areal sweep efficiencies 
at time of breakthrough. At abandonment conditions, the effect of 
mobility ratio on areal sweep is observable, but not to the pronounced 
extent as at breakthrough conditions. 

The conductance ratios show that the majority of the flow 
potential change takes place when the front is near the injection well 
for favorable mobility ratios and near the production well for unfavor- 
able mobility ratios. At time of breakthrough, an apparent change in 
conductance ratio is noted. 

The staggered line drive gave breakthrough sweep efficiencies 
from 4 to 14 per cent higher than those obtained by the direct line drive. 
At abandonment conditions, only small differences existed in areal 
sweep efficiencies for the direct and staggered line drive. However, 
for very poor mobility ratios, over one additional displaceable volume 
of injected fluid was required to reach abandonment conditions with the 
direct line drive. 

The method described for predicting performance of stratified 
reservoirs can be used for any well pattern for which conductance 

ratios and areal sweep efficiencies are available. This data is avail.- 

7 20 

able for the five-spot , the inverted skewed seven- spot , and the 

34 

inverted nine-spot patterns. This study allows similar predictions 
for the direct and staggered line drive patterns. 



5.1 
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The effect of a free gas saturation on areal sweep efficiency 
and conductance ratio has not been sufficiently studied. Accurate 
performance prediction for stratified reservoirs requires additional 
study on the effect of free gas saturations. 



NOMENCLATURE 



a 

c 

d 

Es 

fs 

fw 

h 

k 

L 

i 

M 

Np 

AP 

q 

qw 

rf 

ri 

rp 

rw 

Sew 

Sro 

t 

Vb 

Vd 



well spacing 
constant 

row spacing of wells 

fractional area swept 

fractional flow from swept region 

fractional water flow 

height or thickness 

permeability 

length 

partial length 
mobility ratio 
fractional oil produced 
pressure drop 
total flow 
water flow 
frontal radius 

radius of radial flow region around injection well 

radius of radial flow region around production well 

radius of well bore 

saturation of connate water 

saturation of residual oil 

time 

bulk volume 

displaceable volume - Vb $ (1 - Sew - Sro) 
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v velocity 

Wi water injected 

w width 

a angle open to flow into wellbore from swept region 

y conductance ratio 

$ porosity 

\ mobility - k/gi 

p viscosity 

Subscripts 

A radial flow area around injection well 

B linear flow area between radial flow areas 

C radial flow area around production well 

BT breakthrough 

f field 

i initial 

j any stratum 

m model 

o oil 

s swept region 

t time 

u unswept region 

w 



water 
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.12 .28 .48 .60 .77 1.06 




FIGURE 17. FLOOD FRONT TRACINGS, M = 0. 106 



60 



.09 .29 .47 .67 .91 




.12 .23 .40 .58 .78 1.03 




FIGURE 18. FLOOD FRONT TRACINGS, M = 0.274 
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.09 .24 .52 .90 




Vd 



. 11 



.28 .45 .61 .78 .96 1.27 




Staggered Line Drive 



FIGURE 19. FLOOD FRONT TRACINGS, M = 1 
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.05 .18 .50 





FIGURE 20. FLOOD FRONT TRACINGS, M = 3. 13 
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.11 .25 .40 .60 





FIGURE 21. FLOOD FRONT TRACINGS, M = 12.1 
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TABLES AND FIGURES OF BASIC DATA 
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TABLE III 

DIRECT LINE DRIVE 



Run Initial Pressure Mobility Ratio 

4 1.7 0.106 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 




0.025 


8. 6 


0. 198 


i : o 


- 


0.057 


9.8 


0. 174 


2.0 


- 


0.094 


10. 1 


0.168 


3.0 


- 


0. 135 


10. 2 


0. 166 


5.0 


- 


0.181 


10. 3 


0. 165 


19.27 


- 


0. 672 


11.2 


0.152 


21. 27 


- 


0. 741 


11.6 


0. 146 


24. 27 


- 


0. 809 


12.9 


0.132 


27. 27 




0. 855 


13.6 


0.125 


29. 27 


- 


0. 915 


14.5 


0.117 


31.27 


- 


0.937 


14.9 


0.114 


37.27 


- 


0.954 


15. 3 


0.111 



Note: Areal sweep data is invalid 
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TABLE IV 

DIRECT LINE DRIVE 



Run Initial Pressure Mobility Ratio 

9 2.2 0.106 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


_ 


0. 039 


_ 


_ 


1.00 


- 


0. 104 


10. 4 


0.216 


2. 16 


- 


0. 143 


10.6 


0. 208 


9.75 


- 


0. 340 


10. 8 


0. 204 


13. 50 


- 


0.479 


10.9 


0.202 


20.08 


- 


0. 712 


11.4 


0. 193 


23.41 


- 


0. 844 


11.9 


0. 185 


24.41 


- 


0. 870 


12. 6 


0.175 


25.41 


- 


- 


13.5 


0. 163 


26.41 


- 


0.926 


14. 2 


0. 155 


27.41 


- 


- 


14. 7 


0. 150 


29.08 


- 


0.949 


15. 3 


0. 144 


35.91 


- 


0.974 


15.4 


0. 143 
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FIGURE 22. DIRECT LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 

M = 0. 106 
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TABLE V 

DIRECT LINE DRIVE 



Run Initial Pressure Mobility Ratio 

2 6.5 0.274 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


_ 


0. 107 


10. 7 


0.607 


3.33 


- 


0. 277 


11. 1 


0.586 


6.67 


- 


0.453 


11.5 


0.565 


10.00 


- 


0.697 


12.5 


0.520 


13.25 


- 


0. 860 


14.1 


0.461 


16.58 


- 


0.925 


14.9 


0.436 


23.08 


- 


0.945 


15.0 


0.433 


26.33 


- 


- 


15. 3 


0. 425 


29.67 


- 


0.970 


15. 4 


0.422 



Note: Rerun as run 22 



TABLE VI 

DIRECT LINE DRIVE 
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Run Initial Pressure Mobility Ratio 

9 4,8 0.274 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0. 194 


7. 3 


0.658 


1.0 


1.9 


0. 288 


7.3 


0. 658 


2.0 


3.8 


0. 372 


7.3 


0.658 


3.0 


5. 7 


0.467 


7. 3 


0. 658 


4.0 


7. 6 


0. 558 


7. 3 


0.658 


5.0 


9.5 


0. 665 


7. 3 


0.658 


7.0 


13. 3 


0. 848 


8. 5 


0. 565 


8.0 


15. 2 


- 


9.0 


0.533 


9.0 


17.1 


0.897 


9.2 


0. 522 


10.0 


19.0 


- 


9. 5 


0.505 


11.0 


20.9 


0.930 


- 
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FIGURE 23. DIRECT LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 

M = 0.274 



TABLE VII 
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DIRECT LINE DRIVE 



Run Mobility Ratio 

1 1.0 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 




0.0944 




1.0 


1.0 


- 


0. 240 


C 


1.0 


2.0 


- 


0. 385 


O 


1.0 


3.0 


- 


0. 519 


N 


1.0 


4.0 


- 


0.639 


S 


1.0 


5.0 


- 


0. 726 


T 


1.0 


6.33 


- 


0. 791 


A 


1.0 


7. 33 


- 


0. 840 


N 


1.0 


10. 58 


- 


0.939 


T 


1.0 


17,08 


- 


0.973 




1.0 



Note: Rerun as run 21 
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TABLE VIII 
DIRECT LINE DRIVE 



Run Mobility Ratio 

6 1.0 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0. 271 




1.0 


1.0 


0.84 


0. 356 


C 


1.0 


2.0 


1.68 


0. 450 


O 


1.0 


3.0 


2.52 


0. 546 


N 


1.0 


4.0 


3. 36 


0. 605 


S 


1.0 


5.0 


4. 21 


0. 667 


T 


1.0 


6.0 


5.05 


0. 724 


A 


1.0 


8.0 


6. 72 


0. 790 


N 


1.0 


10.0 


8.42 


0. 845 


T 


1.0 


15,59 


13. 10 


0.934 




1.0 


18. 59 


15.60 


0. 960 




1.0 
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FIGURE 24. DIRECT LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 
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TABLE IX 

DIRECT LINE DRIVE 



Run 

3 



Initial Pressure 
13. 3 



Mobility Ratio 
3. 13 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 




0.099 




_ 


1.0 




0. 210 


- 




2.0 


- 


0. 327 


9. 1 


1.46 


3.0 


- 


0.432 


9. 1 


1. 46 


4.0 




0.482 


8. 1 


1.64 


5.0 


- 


0.530 


7. 3 


1.83 


10.78 


- 


0.654 


7.0 


1.90 


17.40 


- 


0. 750 


7.0 


1.90 


21.25 


- 


0. 798 


7.0 


1.90 


28.00 


- 


0.862 


7.0 


1.90 


34.00 




0. 896 


7.0 


1.90 


35.00 


- 


0.913 


7.0 


1.90 



Note: Areal sweep data is invalid 
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TABLE X 

DIRECT LINE DRIVE 



Run Initial Pressure Mobility Ratio 

8 12.4 3.13 



Elapsed 

Time 


Corrected 

Time 


F ractional 
Area 
Swept 


Pressure 


Conductance 

Ratio 


0 


_ 


0. 051 


_ 


_ 


1.0 


- 


0. 179 


9.4 


1. 32 


2.0 


- 


0. 358 


8. 5 


1. 46 


3.0 




- 


8. 1 


1.53 


4.0 


- 


0. 580 


7.0 


1. 77 


5.0 




0. 642 


6. 5 


1.91 


6.0 


- 


- 


6. 3 


1.97 


7.0 


- 


0. 720 


6. 2 


2.00 


12.08 




0.829 


6. 0 


2.07 


18.58 


- 


0. 920 


6. 0 


2.07 


22.42 


- 


0.952 


6. 0 


2.07 


25.50 


_ 


0.965 


6. 0 


2. 07 
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FIGURE 25. DIRECT LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 
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TABLE XI 

DIRECT LINE DRIVE 



Run Initial Pressure Mobility Ratio 

5 14.7 12.1 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductai 

Ratio 


0 




0.065 






1.0 


- 


0. 101 


- 


- 


2.0 


- 


0. 157 


6.7 


2.19 


3.0 


- 


0. 207 


6. 7 


2. 19 


4.0 


- 


0. 259 


6.7 


2.19 


6.0 


- 


0. 354 


6.7 


2.19 


7.0 


- 


0. 410 


6. 6 


2. 23 


8. 0 


- 


- 


5.9 


2.49 


9.0 


- 


0.496 


5.0 


2.94 


10.0 


- 


- 


4.4 


3. 34 


11.0 


- 


0.525 


3.5 


4. 20 


JL 2.0 


- 


- 


2.9 


5.07 


13.0 


- 


0. 545 


2.8 


5.26 


18. 62 


- 


0.586 


2.7 


5.45 


22.93 


- 


- 


2. 6 


5.65 


29.13 


- 


0. 616 


2.3 


6.40 


39.67 


- 


0.670 


2.3 


6.40 


56.93 


- 


0. 709 


2.3 


6.40 


74.25 


- 


0. 742 


2. 3 


6.40 


91.55 


- 


0. 778 


2. 3 


6.40 


119. 38 




0. 812 


2. 3 


6.40 



Note: Areal sweep data is invalid 



TABLE XII 
DIRECT LINE DRIVE 
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Run Initial Pressure Mobility Ratio 

10 15.0 12.1 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


_ 


0.038 




_ 


1.0 




0. 107 


9.1 


1.65 


2.0 


- 


0. 173 


9.1 


1.65 


3.0 


- 


0. 247 


8.8 


1. 70 


4.0 


- 


0. 323 


8.5 


1. 77 


5.0 


- 


0. 401 


8.3 


1.81 


6.0 


- 


- 


7.0 


2. 14 


7.0 


- 


0. 507 


4.2 


3.57 


8.0 


- 


- 


3.6 


4. 17 


9.0 


- 


0. 551 


3.4 


4.41 


10.0 




- 


3.4 


4.41 


18.0 


- 


0.654 


3.0 


5.00 


24.59 


- 


0. 706 


2.8 


5. 36 


28.34 


- 


- 


2.8 


5. 36 


34.92 


- 


0. 740 


2.8 


5. 36 


41.50 


- 


0. 789 


2. 8 


5.36 


51.92 


- 


0. 807 


2. 8 


5. 36 


62.25 


_ 


0. 833 


2. 8 


5.36 



Run 10 
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FIGURE 26. DIRECT LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 
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TABLE XIII 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

14 2.7 0. 106 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 




0. 122 


11.7 


0. 230 


6.08 


- 


0. 282 


11.8 


0.229 


13.08 


- 


0.482 


12.0 


0. 225 


16.66 


- 


0.598 


12. 1 


0. 223 


23. 16 


- 


0. 768 


12. 5 


0. 216 


26.16 


- 


0.918 


13. 7 


0. 197 


27. 16 


- 


- 


14. 2 


0.190 


28. 16 


- 


0.955 


14. 7 


0. 184 


29. 16 


- 


0.969 


15. 0 


0.180 


30. 16 


- 


- 


15. 2 


0.178 


31.16 


- 


- 


15.4 


0.175 


32.16 


- 


0.978 


15. 6 


0. 173 


37. 16 


- 


0.989 


15. 7 


0. 171 



81 



TABLE XIV 

STAGGERED LINE DRIVE 



Run 

19 



Initial Pressure 

2 . 8 



Mobility Ratio 
0. 106 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0. 0086 


9. 8 


0.284 


1.87 


0. 505 


0.032 


12. 1 


0.231 


5. 20 


2. 21 


0.074 


12.5 


0.224 


12. 12 


5.17 


0.174 


12. 5 


0.224 


18.62 


7.94 


0.279 


12. 5 


0.224 


25. 37 


10.80 


0. 384 


12.5 


0.224 


32.03 


13.65 


0. 503 


12.5 


0.224 


38.62 


16.45 


0. 630 


12. 5 


0.224 


45. 20 


19.20 


0. 717 


12. 5 


0.224 


51.78 


22. 10 


0. 822 


12. 5 


0.224 


55.12 


23. 50 


0. 867 


13. 0 


0.215 


58.45 


24.90 


0.900 


14. 0 


0. 200 


61.70 


26.30 


0.932 


15. 0 


0. 187 


64.95 


27.60 


0. 948 


15.8 


0.177 


74.70 


31.80 


0.975 


16. 3 


0.172 


84.37 


35.90 


0.990 


16. 8 


0.170 
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FIGURE 27. STAGGERED LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 

M = 0. 106 
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TABLE XV 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

12 5.5 0. 274 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


_ 


0.066 


_ 


_ 


2.0 


- 


0. 187 


- 


- 


4.0 


- 


0. 275 


- 


- 


5.0 


- 


0. 347 


- 


- 


6.0 


- 


0. 396 


10. 7 


0.514 


7.0 


- 


0.461 


10. 7 


0.514 


9.0 


- 


0. 584 


10. 7 


0. 514 


11.0 


- 


0. 718 


10. 7 


0. 514 


13. 0 


- 


0.842 


11.3 


0.486 


14.0 


- 


0.923 


11.8 


0.466 


15.0 


- 


0.950 


12.0 


0.458 


17.0 


- 


0.971 


12. 2 


0.451 


21.0 


- 


0.980 


12. 2 


0.451 


25.0 


- 


0.988 


12. 2 


0.451 
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TABLE XVI 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

17 5. 1 0.274 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0.414 


9. 1 


0. 560 


1.0 


0.934 


0. 871 


9.5 


0.537 


2.0 


1.88 


0. 143 


9. 5 


0.537 


3.0 


2. 80 


0.196 


9.6 


0.531 


4.0 


3. 74 


0. 257 


9.6 


0.531 


6.0 


5.60 


0. 372 


9.6 


0.531 


8.0 


7.47 


0.481 


9.6 


0.531 


10.0 


9.34 


0.597 


9.6 


0. 531 


12.0 


11.20 


0. 703 


9.7 


0. 526 


13.0 


12. 10 


- 


9.8 


0.520 


14.0 


13. 10 


0. 793 


10. 1 


0.505 


15.0 


14. 00 


0. 884 


10. 6 


0.481 


16.0 


14.94 


- 


11.0 


0.464 


17.0 


15. 88 


0.951 


11. 3 


0.451 


19.0 


17. 75 


0.976 


11.5 


0.443 


20.0 


18. 68 


- 


11.6 


0.440 


21.0 


19. 60 


0.987 


11. 7 


0.436 



85 




• • • • 
o o o o 

tfl 

W 



<u 



a 

H 



FIGURE 28. STAGGERED LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 

M = 0.274 



TABLE XVII 
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STAGGERED LINE DRIVE 



Run 

11 



Mobility Ratio 

1.0 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0.083 


c 


1.0 


1.0 


1.05 


0. 272 


o 


1.0 


2.0 


2. 10 


0.473 


N 


1.0 


3.0 


3. 15 


0. 685 


s 


1.0 


4.0 


4.20 


0. 823 


T 


1.0 


5.0 


5.25 


0.904 


A 


1.0 


6.0 


6. 30 


0.951 


N 


1.0 


8.0 


8.45 


0.987 


T 


1.0 
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TABLE XVEI 
STAGGERED LINE DRIVE 



Run Mobility Ratio 

16 1.0 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 




0. 113 


C 


1.0 


1.0 


- 


0. 281 


O 


1.0 


2.0 


- 


0.450 


N 


1.0 


3.0 


- 


0.611 


S 


1.0 


4.0 


- 


0. 775 


T 


1.0 


5.0 


- 


0. 876 


A 


1.0 


6.0 


- 


0.940 


N 


1.0 


7.0 


- 


0.967 


S 


1.0 


8.0 




0.992 


T 


1.0 



88 




cn 

w 



Time 

FIGURE 29. STAGGERED LINE DRIVE — AREA SWEPT versus ELA.PSED TIME 
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TABLE XIX 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

13 14.7 3.13 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductanc 

Ratio 


0 


_ 


0.098 


9.1 


1.62 


1.0 


- 


0.229 


8.9 


1.65 


2.0 


- 


0. 369 


8.9 


1. 65 


3.0 


- 


0.509 


8. 8 


1.67 


4.0 


- 


0.636 


8.8 


1.67 


5.0 


- 


0. 739 


8. 1 


1.81 


6.0 


- 


0. 778 


7.3 


2.01 


7.0 


- 


- 


7.1 


2.07 


8.0 


- 


0.832 


7.0 


2. 10 


10.0 


- 


0. 873 


7.0 


2. 10 


13.0 


- 


0.911 


7.0 


2. 10 


15.0 


- 


0.936 


7.0 


2. 10 


17.0 


- 


0.969 


7.0 


2.10 


19.0 


- 


0.988 


7.0 


2. 10 
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TABLE XX 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

18 14.2 3.13 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0. 098 


_ 


— 


1.0 


0. 728 


0.191 


- 


- 


2.0 


1.46 


0. 287 


- 


- 


3.0 


2. 19 


0.377 


- 


- 


4.0 


2.92 


0.465 


- 


- 


5.0 


3.64 


0.553 


10. 5 


1.35 


6.0 


4. 37 


0.632 


10.5 


1.35 


7.0 


5. 10 


0.694 


10.4 


1.37 


8.0 


5. 83 


0. 759 


10. 1 


1.41 


9.0 


6. 56 


0.802 


9.9 


1.43 


10.0 


7. 28 


- 


9.5 


1.49 


11.0 


8.02 


0. 837 


9.2 


1.54 


12.0 


8. 75 


- 


9.1 


1.56 


13.0 


9.48 


0. 866 


9.0 


1.58 


15.0 


10.92 


0.895 


9.0 


1.58 


18.0 


13. 10 


0.920 


9.0 


1.58 


20.0 


14.60 


0.956 


9.0 


1.58 


23.0 


16.80 


1.973 


9.0 


1.58 
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FIGURE 30. STAGGERED LINE DRIVE -- AREA SWEPT versus ELAPSED TIME 
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TABLE XXI 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

15 16.0 12.1 



Elapsed 

Time 


Corrected 

Time 


Fractional 

Area 

Swept 


Pressure 


Conductance 

Ratio 


0 




0.014 


- 


- 


1.0 




0.072 


8.6 


1.86 


2.0 


- 


0. 134 


8.2 


1.95 


3.0 


- 


- 


7.8 


2.05 


4.0 


- 


0. 254 


7. 6 


2. 11 


5.0 


- 


0.315 


7.4 


2.16 


7.0 


- 


0. 424 


7. 1 


2. 25 


8.0 


- 


0. 477 


7.0 


2.29 


9.0 


- 


- 


6. 3 


2.54 


10.0 


- 


0.584 


4.8 


3. 33 


12.0 


- 


0.644 


4.8 


3. 33 


13.0 


- 


- 


4.2 


3.81 


14.0 




- 


3.8 


4. 21 


16.0 


- 


0. 715 


3.4 


4. 71 


18.0 


- 


0.749 


3.4 


4. 71 


24.91 


- 


0. 785 


3.4 


4.71 


42.65 


- 


0. 877 


3.4 


4. 71 


55.08 


_ 


0.915 


3.4 


4.71 
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TABLE XXII 

STAGGERED LINE DRIVE 



Run Initial Pressure Mobility Ratio 

20 16.5 12.1 



Elapsed 

Time 


Corrected 

Time 


F ractional 
Area 
Swept 


Pressure 


Conductance 

Ratio 


0 


0 


0.038 


9.2 


1. 79 


1.0 


0. 523 


0. 104 


8. 2 


2.01 


4. 33 


2. 27 


0. 204 


7.8 


2.11 


7. 58 


3.97 


0.294 


7.8 


2. 11 


10.91 


6. 79 


0.385 


7.8 


2.11 


14. 16 


7.42 


0.468 


7.8 


2. 11 


17.50 


9. 15 


0. 553 


7.8 


2.11 


18. 50 


9. 68 


0.579 


7.8 


2.11 


20. 50 


10. 7 


- 


7.5 


2.20 


21.50 


11. 2 


0. 648 


7.2 


2.29 


22.50 


11. 8 


- 


7.0 


2. 36 


25.00 


13. 1 


0.693 


6.8 


2.43 


28. 25 


14. 8 


- 


6.2 


2. 66 


31.66 


16.6 


0. 747 


5.5 


3.0 


39.25 


20.5 


0. 771 


5.5 


3.0 


48.08 


25. 1 


0.807 


5. 5 


3.0 


57.83 


30. 2 


0. 827 


5.5 


3.0 


67.41 


35. 2 


0. 853 


5.5 


3.0 


76.08 


39.8 


0.872 


5. 5 


3.0 
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FIGURE 31. STAGGERED LINE DRIVE -- AREA SWEPT versus ELAPSED 
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FIGURE 32. DIRECT LINE DRIVE -- MOBILITY RATIOS 

for DISPLACEABLE VOLUMES INJECTED versus 
AREA SWEPT 
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FIGURE 33. DIRECT LINE DRIVE -- MOBILITY RATIOS 

for FRACTIONAL FLOW versus AREA SWEPT 
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FIGURE 34. DIRECT LINE DRIVE -- MOBILITY RATIOS 

for CONDUCTANCE RATIO versus AREA SWEPT 
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FIGURE 35. STAGGERED LINE DRIVE -- MOBILITY RATIOS 
for DISPLACEABLE VOLUME INJECTED versus 
AREA SWEPT 
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FIGURE 36. STAGGERED LINE DRIVE -- MOBILITY RATIOS 
for FRACTIONAL FLOW versus AREA SWEPT 
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FIGURE 37. STAGGERED LINE DRIVE -- MOBILITY RATIOS 
for CONDUCTANCE RATIO versus AREA SWEPT 
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PROGRAM CONRA 

DIMENSION RM( 10) ,DFS(20) ,ES( 10*20) 

READ 1, JE » ( RM ( J ) » J=1,JE) 

READ 1* KE, ( DFS ( K ) » K=1»KE) 

1 FORMAT ( I3»11F7.0) 

DO 10 J= 1 * JE 

2 FORMAT (11F7.0) 

10 READ 2 » ( ESI J *K ) » K=1»KE) 

RM IS THE MOBILITY RATIO AND JE IS THE NUMBER OF MOBILITY RATIOS 
WHICH ARE TO BE STUDIED. DFS IS FRACTIONAL FLOW FROM THE SWEPT REGION 
AND KE IS THE NUMBER OF FRACTIONAL FLOW VALUES TO BE READ IN. 

ES VALUES ARE AREAL SWEEP EFFICIENCIES CORRESPONDING TO THE DFS VALUES. 
SET INITIAL VALUES 
D = 1.0 
A = D 

RE I = 0 «25*D 
RW = D/320.0 

RLW= 0. 5* I D/A - 1.17 + 0 . 63 7*LOGF ( A/RE I ) ) 

ESA = 3.1416*REI**2/(2.0*D*A) 

DO 70 J= 1 * JE 
PRINT 3*RM(J) 

30FORMAT ( 1H1 »40X , 16HMOB I L I TY RATIO =F5.2//10X » 19H CONDUCTIVITY RATI 
10, 1 OX > INFRACTION SWEPT * 13 X , 2HFS/ / 2 3X , 7H1 . 00000 * 1 9X , 1H0 , 1 9X , 1H0 ) 

ESBT = ES { J « 2 ) 



REP = SQRTF( 2.0*D*A*( 1 .O-ESBT ) /3. 1A16 ) 
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DES = 0.01 
ES(J.l) = DES 
DO 70 K=1,KE 

80 IF (ES(J*K) - ESA) 20,20.30 
30 IF (ES(J,K) - ESBT ) 40,50,50 
20 DLW = 0.0 

RF = SQRTF(2.0*D*A*ES(J,K)/3.1416) 

GO TO 60 
40 RF = RE I 

DLW =RLW*(ES( J,K ) - ESA)/(ESBT - ESA) 

DES = 0.05 
GO TO 60 
50 RF = RE I 
DLW = RLW 

600CR = ( LOGF ( RE I /RW ) + 6.28*RLW+ LOGF (REP/RW) ) / ( 1 • /RM ( J ) *LOGF ( RF/RW ) 
1+ LOGF ( REI /RF ) + 6 . 28* ( 1 . /RM ( J ) - 1.0)*DLW + 6.28*RLW+ LOGF ( REP/RW 
2 ) * ( 1 • 0 - DFS(K) )/(1.0-DFS(K)/(RM(J)-DFS(K)*(RM(J)-1.0) ) ) ) 

PRINT 4, CR,ES( J,K) ,DFS(K) 

4 FORMAT ( 10X » 3F20 • 5 ) 

ES(J,K) = ES(J,K) + DES 
IF ( ES ( J , K ) - ESBT) 80,70,70 
70 CONTINUE 



END 



Flow Diagram Program CONRA 





Do Loop 
for all 
M values 







Es bt = Es(j, 2) 

rp = (2da (1 - Es BT ) / TT ) 0,5 

AEs = 0.01 

Es(j , 1) = AEs 







Do Loop \ 
for all ) 
fs values / 
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FORTRAN SYMBOL CONVERSION CHART 



FORTRAN 


N omenc latur e 


Description 


A 


a 


well spacing 


CR 


Y 


conductance ratio 


D 


d 


row spacing of wells 


DES 


AEs 


increment of area swept 


ES 


Es 


fractional area swept 


ESA 


E9 a 


radial flow region around 
injection well 


ESBT 


Es bt 


area swept at breakthrough 


DFS 


fs 


fractional flow from swept region 


DLW 


Jj/w 


partial length to width ratio 


REI 


ri 


radius of radial flow region 
around injection well 


REP 


rp 


radius of radial flow region 
around production well 


RF 


rf 


frontal radius 


RLW 


L/w 


length to width ratio 


RM 


M 


mobility ratio 


RW 


rw 


radius of wellbore 
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FIGURE 38. DIRECT LINE DRIVE -- COMPARISON OF 

ANALYTICAL and MODEL CONDUCTANCE RATIOS 



109 




~ — — Analytical 
Model 



FIGURE 39. STAGGERED LINE DRIVE -- COMPARISON OF 

ANALYTICAL and MODEL CONDUCTANCE RATIOS 
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FIGURE 40 



DIRECT LINE DRIVE -- EFFECT OF L/w RATIO 
ON CONDUCTANCE RATIO 
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EXAMPLE PREDICTION FOR STRATIFIED RESERVOIR 

To illustrate the procedure described for stratified reservoirs, 
assume a staggered line drive unit of 40 acres well spacing composed 
of 2 strata. Under reservoir conditions, the oil viscosity is 4 cp and 
the water viscosity is 0. 4 cp. Assume that no gas saturation exists. 

The reservoir properties are: 





Stratum 1 


Stratum 2 


Height 


20 feet 


20 feet 


Porosity 


o 

« 

LO 

O 


0.25 


Connate water (Sew) 


0. 30 


0. 25 


Residual oil (Sro) 


0. 30 


0. 30 


k Q (at Sew) 


200 md 


75 md 


k (at Sro) 
w ' 


20 md 


22 md 


The bulk volume in each stratum will be the 


same or (20) (80) = 



1600 acre ft. 

Vd = Vb <j> (1-Scw-Sro) 

Vd x = (1600) (0. 30) (1-0. 30-0. 30) = 192 acre-ft 
Vd 2 = (1600) (0. 25) (1-0.25-0. 30) = 180 acre-ft 
The mobility ratios for each stratum are 

M = (V^)/(k o /n 0 ) 

M l = (20/0. 4)/ (200/ 4) = 1 
M 2 = (22/0. 4) / ( 75 / 4) = 3.13 

and 
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k l Vd l _ (200) (192) 

" (75) (180) 



2.84 



To obtain ct for each stratum from equation (25) the inverse 
of the conductance ratio is plotted versus displaceable volumes injected 
(from Figures 13 and 15). The resulting curves (Figure 41) are 
graphically integrated to obtain ct^ and ct^ . (Note that the data for 
each Figure is in the preceding Table. ) 

Cumulative water injected and the corresponding oil produced 
are obtained from Figure 13 in dimensionless terms of displaceable 
volumes injected and fractional area swept. These are then converted 
to Wi and Np and plotted for each stratum versus ct in Figure 42. 
Total water injected and oil produced can be obtained from Figure 42 
by summing the Wi and Np curves for each stratum and is plotted 
in Figure 43. This curve can be converted to an actual time relation- 
ship by knowing the maximum injection rate and well allowables. 
Fractional flow from the swept region is obtained as the slope of the 
Np versus Wi curve and the total conductance ratio is obtained from 
equation (31). Fractional flow is plotted on Figure 43 and conductance 
ratio versus water injected for the system is plotted on Figure 44. 

The above example does not consider the presence of a free gas 
saturation at the start of the water flood. However, it will be within 
the realm of engineering accuracy to assume that oil bank fill up occurs 
simultaneously in all strata. This will be true except in the case of 
high gas saturations because of the unusually favorable mobility ratio 
during oil bank fill up. The volume occupied by free gas prior to oil 
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bank fill up can then be added to the cumulative water injected and 
volume swept to obtain satisfactory performance predictions when there 
is free gas present at the beginning of the flood. 
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TABLE XXIII 



Wi/Vd 

fraction 


Wi 

ac re-ft 


y 2 


!/y 2 


0 


0 


1.0 


1.0 


0. 1 


18.0 


1.48 


0.675 


0. 5 


90.0 


1.50 


0. 675 


0.665 


119. 8 


1. 52 


0.657 


1.0 


180. 0 


1. 75 


0. 571 


1.5 


270. 0 


1.83 


0. 546 


2.0 


360.0 


1. 84 


0. 544 


2.5 


450. 0 


1.85 


0. 540 
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FIGURE 41. INVERSE CONDUCTANCE RATIO versus WATER INJECTED 
for STRATA 1 and 2 



TABLE XXIV 
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Stratum one 


Stratum two 




Wi 


Ct l 


Wi 


( 1 / y) dWi 


ct 2 

2. 84 (1/y) dWi 




0.0 


0 . 0 


0 . 0 


0 . 0 


0 . 0 




100.0 


100.0 


50.0 


35. 38 


100. 5 




200.0 


200,0 


100.0 


68. 83 


195. 5 




300.0 


300.0 


150. 0 


100.88 


285.0 




400. 0 


400. 0 


200. 0 


132. 88 


377.0 








300.0 


188. 18 


535.0 








400.0 


242. 48 


688.0 








500. 0 


296. 38 


843.0 


Wi/Vd Np/Vd 


Wi Np 


Wi/Vd 


Np/Vd 


Wi Np 


fraction fraction 


acre-ft acre-ft 


fraction 


fraction 


acre-ft acre-ft 


0. 769 


0.769 


147. 5 147.5 


0.665 


0.665 


1 19; 5 119.5 


0.80 


0. 795 


153.5 152.5 


0.80 


0. 743 


143. 8 133.5 


1.0 


0. 90 


192.0 173.0 


1.0 


0.803 


180.0 144.2 


1.2 


0.952 


230.0 182.8 


1.5 


0. 883 


269.0 158.7 


1.4 


0.985 


269.0 189.0 


2,0 


0.930 


359,0 167.0 








2.5 


0.961 


449.0 173.0 
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URE 42. WATER INJECTED and OIL PRODUCED versus ct 
for STRATA 1 and 2 
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TABLE XXV 



Wi x 


wi 2 


Np x 


Np 2 


Wi 


Np 


100 


50 


100 


50 


150 


150 


145 


75 


145 


75 


220 


220 


200 


105 


173 


105 


305 


278 


225 


120 


180 


120 


345 


300 


250 


132 


185 


127 


382 


312 


300 


160 


190 


138 


460 


328 


350 


187 


190 


146 


537 


336 


400 


216 


190 


152 


616 


342 


500 


276 


190 


160 


776 


350 


600 


344 


180 


166 


944 


356 


700 


410 


190 


171 


1110 


361 



Wi 


fo 


fw 


220 


1.000 


0.000 


300 


0.590 


0.410 


340 


0. 550 


0.450 


350 


0.400 


0.600 


400 


0.245 


0. 755 


500 


0. 127 


0. 873 


600 


0,054 


0.946 


800 


0.038 


0.962 


1000 


0.031 


0.969 



500 
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FIGURE 43. OIL PRODUCED and WATER CUT versus WATER INJECTED 
for STRATIFIED SYSTEM 
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TABLE XXVI 




k 2 Vd 2 



k 2 Vd 2 
+ ^T 



^2 



1 + 0. 352 y_ 

r w 



Wi 2 


!/y 2 


Yl 


Wi 


Y 


25 


0. 675 


1.480 


75 


1.126 


50 


0. 672 


1.490 


150 


1.128 


78 


0. 670 


1.492 


225 


1. 128 


105 


0. 660 


1. 518 


305 


1. 134 


160 


0. 600 


1.669 


460 


1.172 


216 


0. 552 


1. 812 


616 


1. 210 


276 


0. 550 


1. 820 


776 


1. 214 


344 


0. 542 


1. 845 


944 


1.220 


410 


0. 540 


1.850 


1110 


1. 221 
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FIGURE 44. CONDUCTANCE RATIO versus WATER INJECTED 
for STRATIFIED SYSTEM 
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